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Abstract 

The total cross sections of the pp — > K + Ap reaction near threshold 
measured at the COSY(Juelich) synchrotron using the internal target 
facility COSY-11 are described in the frame of the model taking into 
account the one pion and one kaon exchange as well as the final state 
interaction mechanisms. It is shown that near the threshold all these 
mechanisms give the compatible contributions, but with the increasing 
energy the one kaon exchange mechanism becomes the dominating 
one. 

Introduction 

The total cross sections of the pp —* K + Ap reaction near threshold measured 
at the COSY(Juelich) synchrotron using the internal target facility COSY-11 
[1, 2, 3] present good test for the various models of the associated strangeness 
production. In Tab. 1 the 13 values of cr to t(pp — > KAp) measured at different 
periods [1, 3] are presented. Let us dwell on some theoretical models of the 
pp — ► K + Ap reaction. In the paper [4] the parameters of the model presenting 
the amplitude coherent sum 

M = Mqpe + Mqke + M FS i {7t) + M FS i {K ) , (1) 
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Table 1: Total cross sections a tot {pp — > pKA) from [1] and [3]. Q = W — 
m p — m A — m K is the excess energy. 

had been chosen in order to match the existing total cross sections in the 
energy interval 2 < T p < 6 GeV [6] 1 . It had been shown, that the main 
contribution to the total cross sections in this energy interval provides the 
one kaon exchange mechanism (OKE). Final state interaction (FSI) con- 
tribute very little except the region close to the kaon production threshold. 
The Juelich group [7] also has stated that, in principle, K exchange alone 
could explain the total cross section especially after inclusion of FSI effects. 
However the 7r° exchange cannot be neglected near the threshold. It was 
argued that the experimental data require a destructive interference between 
7T° and K exchange contributions. In the paper [8] the detailed study of the 
chances for identifying the reaction mechanism of strangeness production are 
performed by considering the 7r° and K exchange diagrams. 

The aim of present work is to perform the calculations as accurate as 
possible using as the amplitude a coherent sum (1), thus following the papers 
of J.M.Laget [4, 5]. The model, which describes fairly well the total cross 
section in the wide interval of the excess energy Q = W — m p — m\ — m K 
from 0.7 to 500 MeV, is planned to use for the analysis of the proton and 

1 In the next Laget's paper [5] the data close to threshold, (0.68 < Q < 6.68 MeV [1]) 
are already taken into account. 
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kaon spectra in pp — > pAfT reaction at COSY energies [9]. 

Although near the threshold the spin flip effects are negligible, the special 
attention is paid to the complete taking into account of the spin dependence 
of the considered amplitudes. It allows to make realistic predictions of the 
polarization observables for intermediate energies, which will be published 
elsewhere. The Stapp formalism of the M-functions is the most proper for- 
malism in this case and it is shortly described in the Appendix A. 



1 The total cross section 

The five-fold cross section is equal to 

dV - Pfe 2 * I I | M | 2 , (2) 

dppdVLpdVLx (2tt) 5 8E p Ap b m \p K E A -p A E K z AK \ 

where z A k = Pk-Pa/pkPa- The amplitude M has been also calculated as a 
coherent sum (1) where the separate contributions in this sum are considered 
in details in the next sections 

From eq.(2) the total cross section reads as 

atot = 1 1 / d V / dn K ^- -\M\ 2 . (3) 

(2tt) 5 32mp b J Pp J K E P \p K E A - p A E K z AK \ 11 1 ; 

Let pb be the momentum of the beam proton, pb being directed along the 
z-axis. The total 4-momentum and the squared invariant mass of the system 
are equal 

p tot = (E, p 5 ) = (m p + ^m 2 p + p 2 b , p b ) , s = W 2 = p 2 tot . 
The proton momenta are placed then inside the ellipsoid 



P 2 px + Ply (Ppz-Po) 2 



where 



2 (s - (m p + m A + m K ) 2 )(s - (m p - m A - m K ) 2 ) 

is squared cm. momentum of proton in the case of minimal mass of KA 
sub-system equal to m A + uik and 

s + ml - (m A + m K f 
Po = Pb — • (6) 
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It is natural to prepare the integral (3) for numerical integration over d 3 p p 
by choosing the elliptical coordinates 

p px = a sinh u sin 9' cos <f> 

P P y = a sinh u sin 9' sin <p (7) 
Ppy — Po + a cos h u cos 9' , 



where 

and 
with 



_ Ptot 

W ' 



0<u< u max 



sinh u m ax — ) cosh u max — 

Ptot Ptot 



The Jacobian of transition from Cartesian to elliptical coordinates is equal 
to 

d 3 p p = a 3 (sinh 2 u + sin 2 9') sinh u sin 9' dwd^'d^ 



or 

3/_ i2„, 12 



d p p = a (cosh u — x )d coshwdx'd^ . (8) 
Note that the volume of this ellipsoid is equal to 

t t- 4tt 3 E 
3 Hcm W 

Indeed it is calculated with the Jacobian (8): 

d 3 p = a 3 / dcosh-u dx d0(cosh 2 u — x 2 ) = 

E 

o I Ptot 9 2. All n E 

2na 3 d cosh M (2 cosh 2 - -) = -j&™w 

We replace the measure dflx, where the angles Q = (9k, 4>k) of kaon are 
defined with respect to fixed frame, by dflK(KA) with the kaon polar angles 
defined with respect to momentum of KA sub-system. In this case we can 
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make the trivial integration over azimuthal angle of proton and obtain from 
eq.(3) 



I l rE/ptot rl r 

atot = ^32^/1 ^ujjx'j dn K{m 

r? I Ml 2 

(cosh 2 .-^ m r . (9) 

£> P \Pk&a — Pa&kZak I 

Let the variables of integration in eq.(9) be replaced by variables i = 
1, 2, 3, 4 such that ^ qi ^ 1. Then 

cosh-u = 1 + ( l)gi , dcosh-u = (— — l)dgi , 

Pb Pb 

x' = 1 + 2q 2 , da:' = 2dg 2 , 

x K = 1 + 2q 3 , drr^- = 2d<? 3 , 

0x = 27rg4 , d(px = 2ndq4 . 

Then we can rewrite the eq.(9) as 

<?tot = , Q \, Q a (— - 1) / dgi / dq 2 [ dq 3 [ dq A 

(cosh 2 u - x' 2 )^ , - lMl \ r . (10) 

E p \p K E A - p A E K z AK \ 

With the eq.(10) we have tested the model considered in next sections 
by using the CERN program RIWIAD, an adaptive multidimensional inte- 
gration subroutine which permits numerical integration of a large class of 
functions, in particular those that are irregular at the border of the integra- 
tion region. 



2 Pion exchange mechanisms 

In what follows the tensor-like notation for the amplitudes with the indices 
characterizing spin projections is used. The upper and lower indices relate 
to the final and initial channels, respectively. The same index at the up and 
down positions stands for a summation over this index. So, the pp — * pAK + 
amplitude looks like 

M%Z b A ( P K,PA,P P ;Pt,Pb) , 
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Figure 1: OPE (a,b) and FSI diagrams (c,d) 



where at, a b , a p and a A are the spin projections of the target, beam, final 
proton and A respectively. PK,P\,P P ,Pt,Pb are the corresponding momenta. 



2.1 One Pion Exchange amplitude 

The one pion exchange(OPE) graphs are shown in Figs. 1 (a,b). Let us start 
with the contribution of the graph a: 

M {7T o p ^ K+A) a ^(p A ,p K ;p t ,p b - P P )G {n o pp) l P b (p p ,p b - p p ;p b ) 

t — m* 

where t = (p b — p p ) 2 < is the squared mass of the virtual n°, M^o p ^k+k) 
is the half-off-shell 7r°p K + A amplitude and G is the n pp vertex which is 
equal in Stapp formalism to (see Appendix A) 

G^ PP C b ( Pp ,Pb-p P ;Pb) = ^f(t)m(e - V p V b )H , (12) 

where g n is the ttNN coupling constant, and f(t) is the pion form factor for 
which we take the monopole representation 

m 2 _ A 2 

/(*) = » ( 13 ) 

A,r being the cut-off parameter. Defining the pion 'wave function' in the 
nucleon 

a v _ G (^vp)7 h = mg^ ( 1 1 

t-ml ~ y/2\t-ml t-Al 
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and replacing the amplitude M^o p _^ K + A ) by the explored in the Rutherford 
Laboratory amplitude Mr w - P ^K A) [10] , which are connected as follows 

we rewrite the eq.(ll) as 

^M {n - p ^ K o A f^(p A ,p K ;p t ,p b - p p )^ {n o pp) l P b (p p ,p b - p p -,p b ) . (15) 

The contribution of the graph b pictured in Fig. 2 is obtained from the 
expression (15) by the interchange of the beam and target quantities and we 
derive the final expression of the OPE contribution 

M {0 PE) (Pp ,PA,PK]Pb,Pt) = 
-j=M {7 r- p ^ K o A) ^(p A ,p K ;p t ,p b -Pp)^{nOppf a P b ipp,Pb -PplPb) ~ 

(&<-*), (16) 

where the relative minus sign is required by Pauli principle. In the Appendix 
A it is explained how the M-function expression (12) of the 7r°pp vertex is 
derived and how to obtain the M-function of the reaction n~p — > K°A from 
the canonical flip and non-flip complex amplitudes G s and H 8 . We have 
taken into account only the S- and P-wave partial amplitudes in the partial 
wave development [11] 

1 oo 

1 oo 

H' = -Y,[T»-Ti-\Pfo)i ( 17 ) 
9 i=i 

where q = y^qj, qij are the cm. momenta in initial and final channels 
and x is the cosine of the cm. scattering angle. Following the papers 
[10](Bakerl978) we took 



Table 2: Resonance parameters of the n p — > K°A amplitude 
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where 



21+1 



2m r 1 + (1.7q r ] 



V^ + m r l + (1.7g) 2 (18) 

and the parameters are given in Tab. 2. Unfortunately the complex coeffi- 
cient a describing the background is not presented by Rutherford Laboratory 
group, so we were to fit it to describe the differential cross sections at four 
cm. energies 1.661, 1.683, 1.694 and 1.724 GeV given in [10] (Baker 1978). 

The virtuallity of the initial n° is taken into account as follows. At first 
the amplitudes (17) are multiplied by the same form factor (13) as the ir°pp 
vertex (12) 



G" 




J2i(i + w+ + m-]Pi(x) , 



£ [T l+ - T,_] Pfts) 



(19) 



In addition the cm. momentum in initial channel 7T°p is calculated with the 
virtual mass of the pion t — (p& — p p ) 2 < 0: 



Qi 



V[(Wka ~ m p y - t][(W KA + m p y - t 



2W KA 

where Wka is the invariant mass of the reaction. 



(20) 



2.2 FSI mechanism with the intermediate pion 

The FSI graphs are shown in Figs. 1 c, d. Writing the contribution of, for 
example, the graph c in accordance with Feynman rules and putting the 
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proton 'spectator' p s on the mass-shell we obtain 



2 



M FS i ilt)t 3£ = -i j j^^ s 



1 

m p f dp s 



M (paAv ^ pA) ^— 2 M { ^p t ^K+K)l v t ®{^p P )l s b 



]p_ f c 

pA J (2tt 



W pA J (2Ttf2E s 

M {PaKv ^ pK) 2 _ M { ^ Pt _ K+Kv) ll $ ( ,rV)% , (21) 

where is the mass of virtual A inside the triangle loop, the pion 'wave 
function' in the nucleon $( n o pp ) is defined in eq. (14) and the non-relativistic 
approximation of the A propagator 

m 2 2 ^ ^pA / 2 _^2n 
m v m A —Won S J 

/ / L 

is used. Here W^,a = y/Sp~A, s p a = p 2 pA an d 



\/[(W p a - m p ) 2 - m A ][(W P A + m p ) 2 - m\ 



2W pA 

^ = Vi(W P A - m p ) 2 - ml] [{W pk + m p f - ml] 

2W p a 

are the cm. momenta of A after (on-shell) and before (off-shell) the scatter- 
ing. 

Simplifications are necessary in order to calculate the integral (21) since 
it requires a knowledge of the off-shell ir°p — > K + A and pA amplitudes. In 
this we have followed the paper [12] with some modification described in 
[13]. Namely, these amplitudes are taken out of the integral sign at some 
momentum p° s placed on the singular surface of the integral corresponding to 
the mass-shell of the virtual A. 

When a separable potentials for pA scattering are assumed for each partial- 
wave state the relation 

M° p { f = f{ml)M;i (23) 



2 As it was explained above, the same index a v , spin projection of the virtual A, at the 
up and down positions stands for a summation over this index. 



9 





Set A-l 


Set A-2 


Set A-3 


Set B-l 


Set B-2 


'So 


0.209 


0.256 


0.239 


0.257 


0.258 




0.217 


0.246 


0.252 


0.247 


0.243 



Table 3: Parameter /3(GeV/c) of separable pA potential by Y.Takahashi et 



is strictly valid. The usual form of the form factor / is 



f(m 2 v ) 



q 2 

T-c.m. 



+ P 2 



e+P 2 



(24) 



where (3 is the cut-off parameter. Near threshold two S-waves, 1 Sq and 3 S'i, 
give the major contribution. In the paper[14] different sets of separable pA 
potential were found with good fits to the cross- sect ion. They are presented 
in Tab. 3. Following this table we neglected the difference between the spin 
singlet and spin triplet S-wave parameters and took the form factor (24) for 
the whole low-energy amplitude with f3 varying in the region 0.1 < (3 < 0.3 
GeV/c. Taking into account that 
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the eq. (21) can be rewritten as 
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M W ~ i mp f dPs 



1 1 

+ 



(25) 



We took the half-off-shell 7c°p — > X+A and on-shell part of pA amplitudes 
(23) out of the integral sign at specially chosen momenta. Following the ar- 
gumentation of [13] the best choice for the graph c is the maximal momentum 
of the on-shell loop proton p s , which is directed along p A and is equal to 

F! cm TP „cm 

p p(+ ) = pm(^ + ^— )• 

In this case A is on-shell too and has correspondingly the minimal momentum 

PA(-) 

Tpcm TP „cm 

/ Aa_ -fcpA q x 

Pa(-) = PpaIttt T7T- —) ■ 

VVpA WpA PpA 

Here E^ m and E A m are cm. energies of the proton and A respectively. Under 
these assumptions we can rewrite eq.(25) as follows 
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M FS I(n)cl^(Pp,PA,Pk;Pb,Pt) = 

-iM {pK fj2{p P ,PK\P P {+),PK{-)){e - V p(+) V b )l 
M {7T o p ^ K+A) l v t (p K ,p A{ _y,p b -p p{+) ,p t ) ■ F(s pA ,t) , 

where t — (pt — Pk) 2 and 



(26) 



F(s pA ,t) 



W pA V2j (2nf2E s 



(Pb-Ps) 2 - ml {p b ~p s ) 2 -Al 



+ 



(27) 



In Appendix B the details of the calculation of the integral (27) are given. 
The result is 



3 

£j_ + a'j aj aj ' ' 2((3 + aj)\ ' 

where q b is the cm. momentum of the system p b + (p t — p K ) equal to 



—% In 



t 



3+ -r- a j 



32nW pAqb V2 

— 2 arctan — 2 arctan + 4 arctan 



(28) 



y^Pg + m) 2 ~ t}[(W pA - my~F] 
2W pA 



t = (pt~ Pk)' 



(29) 



and 



mz 



{ 1± = *(i-^)± 9 



cm. ; 

Ci = l, 



«2 = V 1 " " ' 

G± = — ^2) ± 9c.m. 5 

C 2 = -1 . 



(30) 
(31) 

The contribution of the second FSI graph with virtual pion (Fig. 1 d) is 
derived by the interchange b <-> t: 

M FS I(n)d%%(Pp,PA,Pk;Pb,Pt) = 

-iM (pA) ^(p p ,p A ;p p( _ ) ,p A(+) )(e - V$, ( -)^)£ 
M {n0p ^ K + A) l v t (p K ,p A{+) ;p t - Pp^^Pb) ■ F(s pA ,t) , (32) 
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Figure 2: OKE (a,b) and FSI diagrams (c,d) 



where t = (j>k — Pb) 2 , and 

fpcm 7-1 cm rpcm rp „cm 

Pp(-) = PMI777 TFT ~ / ' Pa(+) = Pm(ttt 77T~ —J ■ 

l^pA l^pAPpA l^pA l^pA PpA 



3 Kaon exchange mechanisms 
3.1 One Kaon Exchange amplitude 

The OKE graphs are shown in Figures 2a,b. The final expression of the OKE 
contribution is similar to eq.(16) 
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M {OK E)l A b ^ (Pp ,PA,PK]Pb,Pt) = 

iM {K+p f^(p p ,p K ;p t ,p b - p A )$ {K+pA f a A (p A ,p b - p A ;p b ) - (b <-> t) . (34) 
The kaon 'wave function' in the K + pA vertex is equal to 

*W)* A = qkV^ {t^f k - t=-£) (e " ^ aH) " ' (35) 

where t = (p A — Pb) 2 , 9k is the K + pA coupling constant [15] and Ak is the 
cut-off parameter in the monopole form factor (see eq. (13) for virtual pion) 
of virtual K. 

The M(x+ P ) is the half-off-shell K + p K + p amplitude, the M-function of 
which is derived from the canonical amplitude in the same way as M^- p ^ K o A ) 
(see the previous section and Appendix A). In order to calculate the canon- 
ical flip and non-flip amplitudes we have utilized the isotriplet partial-wave 
amplitudes below 1.5 GeV/c obtained in the phase shift analysis by Martin 
[16]. 

3.2 FSI with intermediate kaon 

For the graph in Fig. 2 c we have 

M F SI{K)cl A b ° P t {Pp,PA,Pk]Pb,Pt) = 

- iM (M)"l^(P P ^A;p P (-),PA(+)) (e - V A{+) V b )ll 
M (K+ P C t (PK,Pp(-y,Pb-PA(+),Pt) ■ F K (s pA ,t) , (36) 

where t = {j>k — Pt) 2 and 
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F(s pA , t) = 



g K ni A y/mm A 
W pA 



dp. 



(2tt) 3 2£ s 



(Pb - Ps) 2 - m 2 K (p b - p s ) 2 - A 2 K 



1 



+ 



(37) 



where now the integration is performed over the 3-momentum p s of the on- 
shell A in the triangle loop. Manipulations similar to those described in 
Appendix B for the pion case result in 
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2(0 + a,) 



(38) 



where 



ai = E b \ 1 



(m 2 + m\ — m 2 K ) 2 



4m 2 m A 



«2 = E b \ 1 



(m 2 + m\ — A 2 K ) 2 



i.m 2 m\ 



m 2 + m 2 A -m 2 K 

?i± — lb ~ ± q c 

2mm A 



m 2 + m 2 A -A 2 K 
?2± — % ± 1, 



(39) 



2mmA 

Ci = 1 , C 2 = -1 , (40) 



and q 



q b . The contribution of the second FSI graph with virtual kaon 



(Fig. 2 d) is derived by the interchange b <-> £. 



4 Results and discussion 

Let us start with the total cross sections for OKE+FSI mechanism. Free 
parameters in this case are gx - K + pA coupling constant, Ax - the cut-off 
parameter in the monopole form factor of virtual K and j3 - the cut-off param- 
eter in the form factor in half-off shell pA interaction. The cut-off parameter 
Ax = 0.9 had been chosen and with g K = -4.l7y/4n = -14.78 GeV, the 
value recommended in the paper [15], it appeared that the experimental data 
presented in Tab. 1 are fairly well described by only OKE mechanism. Still 
to find place for the FSI contributions, from the one hand the coupling con- 
stant gx had been minimally decreased to the value gx = —11.8 in order to 
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stay in agreement with SU(3) predictions for the kaon-hyperon-nucleon cou- 
pling constants (see [15]), and from the other hand the FSI contributions had 
been minimized by choosing rather small value of the half-off shell pA cut-off 
parameter: (3 = 0.1. The resulting total cross sections versus the excess en- 
ergy Q = W — m p — itia — itlk are shown as double logarithmic plot in Fig. 



3. It is seen that OKE+FSI(K) model alone with the plausible parameters 



pp — > pKA), mb 




10 



10 



g pKA =-11.8, A K =0.9GeV 

j i i i l_ 



10 2 10 1 

Q=W-W lh , GeV 



Figure 3: Red circles are the experimental total pp — > pAK + cross sections 
versus the excess energy Q = W — W t h- Green curve presents the calculations 
in the frame of OKE model. Blue curve - the calculations in the frame of 
OKE+FSI model. 

describes very well the experimental data, which is already a hint in favor 
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of kaon mechanism dominance in pp — > pAK + reaction. The Q behavior of 
the difference between pure OKE and OKE+FSf(K) curves demonstrates, 
that near the threshold FSI contribution is compatible with the OKE one, 
whereas at higher Q the OKE becomes the dominating mechanism. 

Though the kaon mechanism seems to be enough we should find place 
for the OPE+FSI(7r) graphs. Free parameters in this case are g n - the nNN 
coupling constant and A„. - the cut-off parameter in the monopole form factor 
(13) of virtual n. The cut-off parameter in the form factor in half-off shell 
pA interaction (3 = 0.1 was fixed by OKE+FSI(K) fit. The well established 
value g n ~ 13.5 (g^/4ir = 14.5) had been fixed. Then it appeared that using 
the common accepted A ff ~ 1 GeV results in high OPE contribution near the 
threshold which would be impossible to compensate by other mechanisms. 
We have decreased the cut-off parameter A T to the value A w = 0.515 GeV. 
The resulting plots are shown in Fig. 4. It is seen that only for Q < 1.7 MeV 
the contribution of OPE+FSI(7r) graphs exceeds experimental data and there 
is no possibility to decrease this contribution. For higher Q the OPE+FSI(7r) 
contribution to the experimental data is negligible. 

The applying the full model, including the both kaon and pion mecha- 
nisms, with the parameters fixed above had shown that in case of " construc- 
tive" interference 

M = Mqpe + M OKE + M FSI(jl) + M FSI{K) 

the model total cross sections exceed the measured ones, the difference being 
especially large near the threshold. In contrast, the choice of destructive 
interference 

M = -Mope + M OKE - M F5 / W + M FSI{K ) 

results in fairly well description of the data, which is demonstrated in Fig. 
5. This our observation confirms the conclusion of the Juelich group [7] that 
the experimental data require a destructive interference between n° and K 
exchange contributions. 

In conclusion we can confirm the main features of the pp — > pAK + scatter- 
ing at low and intermediate energies observed in the previous studies [4, 7, 8] : 
the kaon mechanism dominance at intermediate energies, the compatibility of 
all mechanisms near the threshold and the destructive interference between 
kaon and pion amplitudes. The model presented in this work describes suffi- 
ciently well the measured total cross sections pp — > pAK + in the wide interval 
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of the excess energy Q = W — m — m\ — mx- 0.68 < Q < 430 MeV. Our 
further plans are to apply this model for the description for the analysis of 
the proton and kaon spectra in pp — > pAK reaction at COSY energies [9] 
and to make predictions for the polarization observables at these energies. 
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A M-functions 

When calculating a coherent sums of various mechanisms one needs in prin- 
ciple to transform the matrices of corresponding amplitudes from center of 
mass frames to the common laboratory frame. These transformations depend 
on the spin basis chosen for the one-fermion states. Usually the canonical 
or helicity basis are used which are transformed with unitary two by two 
matrices depending on a fermion momentum. They are so called Wigner 
rotations. The covariant basis [17]transforming independently of the par- 
ticle momentum by the unimodular two by two matrices is free from this 
deficiency. Amplitudes in the covariant basis are usually referred to as the 
M-functions of Stapp. The covariant formalism of the M-functions developed 
by Stapp [17] is based on the matrices <r M = (1,<?) and cx M = (1, — <r), where 
a are the standard Pauli matrices. With each momentum P M of the reaction 
two by two matrices P = P^o^ and P = P^o^ are associated. The products 
PiPj of these matrices are the elements from which the M-functions are built. 
The matrices V* and Vi associated with the four-velocity of the i-th particle, 
serve as the metric tensors of the particle when performing the contraction 
over this index (traces, successive processes). 

Let us explain the M-function expression (12) of the ir°pp vertex and how 
to obtain the M-function of the reaction n~p — > K°A from the canonical 
amplitude. As to the n°pp vertex it is nothing but the usual expression 
gUf^Ui with the exception that the modified Dirac bispinors are used. In 
the Weyl representation they are equal to 




yba ^ 



(41) 
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In the Weyl representation the 75- matrix looks like 

/ e 
75= ^0 -e 

and the amplitude of transition between real states with the production of 
the pseudoscalar is equal to 

K(Pf, 9; P l ) = 9Kl5 pu p b = 9y^Jm(e a b - VfV m ) , 

or in indexless form 

M(p f ,q;pi) = gu f ~f 5 Ui = g^m f mi(e - V f Vi) , (42) 

which explains the eq.(12). 

As to the M-function of the 0| — > 0| reaction it is possible to represent 

as 

M = G m h + H m b 2 , (43) 

where G m and H m are the complex functions and b 1>2 are the basis M- 
functions equal to 

b = e + VfVj b ^ a(e + gjg + (V f V + VVj) 
1 a/ (Mi, Mi) ' 2 V / (M 2 ,M 2 ) 

Here V, V* and V) are the four- velocities of the cm., of initial and of fi- 
nal spinor particle. The coefficient a and the normalization coefficients 
a/ (Mi, M^ are equal to 

o (V.*S+Vf) 

(m 1 ,m 1 )= 2[i + (y / ,y i )] , 

/ T\/r -\ir \ o [i-(v.V/) 2 ][i-(v.v0 2 ]-[(V/.Vi)-(v.V/)(v.Vi)] 2 

(M 2 , M 2 J - ^ l+^.VO • 

In this case the corresponding cross section is equal to |G m | 2 + |iJ m | 2 . Thus 
there is the unitary connection between G m , H m and the canonical flip and 
non-flip cm. S- matrix amplitudes 

S = G' + iH"h, h = n k e k , n~ h * kf 



\ki x kf\ 
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where kij are initial and final cm. momenta [11]. It is 



G m = cos G s - sin H s , 
2 2 

H rn = sin G s + cos £ # s , 
2 2 



where the angle f is determined by 



j, I (Jo — (V,V,) + \{V,Vi) + \ „ 

e ' = fe^' -* s V (^/)-Kv,vt)-i - " e • (45) 

The 6* is the cm. scattering angle. 



B Calculation of the transitive form factor F 

Let us chose the pA cm. system. Then 



p b = (E b , q 6 ) , p s = f , (p 6 



2m 2 - 2£ fe £ 5 + 2g 6 ^a: 



E h 



s p a + m 2 -t 



2W, 



<lb 



vP> + - *] PV - ™) 2 - *] 



pA 



2W 7 , 



pA 



Choosing z-axis along q&, taking into account d£ = d0d:r£ 2 d£ and integrating 
with respect to azimuthal angle and x we derive 



F(s pA ,t) 



9tT 



y/216ir 2 q Jo 



In 



2E h E i - 2q b £ - 2m 



1 



+ 



1 



ml 



In 



2E b E ii - 2q£ - 2m 2 + A 2 



2E h Ei. + 2q h i - 2m 2 + ml 2E b E i + 2q£ - 2m 2 + A 2 



(46) 



where q = ^^qb- This integral can be computed if we replace the arguments 
of the logarithms by the parabola with the roots coinciding with the roots of 
these arguments. They are complex and equal to 



£(m 7r )± = q b (l 



2m 2 



±iE, 



A 2 

1 - (1 ^) 2 

v 2m 2 ' 



(47) 
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and we have to calculate 

9n 



F(s p a, t) 



647r 2 g\/2 



1 



1 



+ 



1 



+ 



£ - q cm - it £ + q C m + it £ + i/3 £ - i(3 



ln e-e(A,) + ln e-e(A.) 



Here arguments of the each logarithms has singularities only in the lower or 
upper half-plane. The first and the third logarithms are free from singularities 
in lower half-plane and closing the contour for them in lower half-plane we 
meet the poles at — q cm — ie and — ij3. In contrast the second and forth 
logarithms pick up in upper half-plane the poles at q cm + ie and ij3. This 
yields 



F{s pA ,t) 



9n 



(-In 



-qcm-ic + timn)- 



^ -q cm -ie-£{K) + 
-q cm -ie + ^(A 7r )_ 



9-k ; 
167rg 



-iP + ^rriir). 

ln-^'^ + +ln 

-iqcm + ikVn^)^ 



-jqcm - i£(Air)- 

-iqcm + itiAn)- 



(49) 



+ 



13 + i£(m„) 



Defining 



a.1 = Ei 



m 



IL)2 



i± 



?6(1- 



2m 2 



2m 2 ' 



a 2 = E h \l 1 

6± 



^2 

2m 2 / 



- ^2 ) ± g c 
C 2 = -1 , 



Ci = l, 

we can write the final expression in the form similar to that of Laget 



(50) 
(51) 



F(s pA ,t) = 



—% ln 



e ? 2 _ + « 2 



— 2 arctan 



<5j+ 
a, 



— 2 arctan 



6"- 



ft.; 



+ 4 arctan 



2(0 + a ) 



(52) 
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cr lot (pp —> pKA), mb 




Q=W-W th5 GeV 



Figure 4: Red circles are the experimental total pp — > pAK + cross sections 
versus the excess energy Q = W — W t h- Green curve presents the calculations 
in the frame of OPE model. Blue curve - the calculations in the frame of 
OPE+FSI(vr) model. 
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Figure 5: Red circles are the experimental total pp —> pAK + cross sections 
versus the excess energy Q = W — W t h- Blue curve - the calculations in the 
frame of full OKE-OPE model (destructive interference) 



25 



